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ABSTRACT 

We use the abundance and weak lensing mass measurements of the SDSS maxBCG cluster catalog 
to simultaneously constrain cosmology and the richness-mass relation of the clusters. Assuming a flat 
ACDM cosmology, we find cr8(^m/0.25) 041 = 0.832 ±0.033 after marginalization over all systematics. 
In common with previous studies, our error budget is dominated by systematic uncertainties, the 
primary two being the absolute mass scale of the weak lensing masses of the maxBCG clusters, and 
uncertainty in the scatter of the richness mass relation. Our constraints are fully consistent with the 
WMAP five-year data, and in a joint analysis we find a s = 0.807 ± 0.020 and Cl m = 0.265 ± 0.016, 
an improvement of nearly a factor of two relative to WMAP5 alone. Our results are also in excellent 
agreement with and comparable in precision to the latest cosmological constraints from X-ray cluster 
abundances. The remarkable consistency among these results demonstrates that cluster abundance 
constraints are not only tight but also robust, and highlight the power of optically-selected cluster 
samples to produce precision constraints on cosmological parameters. 

Subject headings: cosmology: observation — cosmological parameters — galaxies: clusters — galaxies: 
halos 



1. INTRODUCTION 

The abundance of galaxy clusters has long been rec- 
ognized as a powerful tool for constraining cosmological 
parameters. More specifically, from theoretical consid- 
erations (e.g. Press & Schechter 1974; Bond et al. 1991; 
White et al. 1993; Sheth & Tormen 2002) one expects the 
abundance of massive halos to be exponentially sensitive 
to the amplitude of matter fluctuations. Though some 
theoretical challenges remain (see e.g. Robertson et al. 
2008; Stanck et al. 2009), this basic theoretical predic- 
tion has been confirmed many times in detailed numerical 
simulations, and a careful calibration of the abundance 
of halos as a function of mass for various cosmologies 
has been performed (see e.g. Jenkins et al. 2001; Warren 

Center for Cosmology and Astro-Particle Physics, The Ohio 
State University, Columbus, OH 43210 erozoOmps.ohio-state.edu 

2 Kavli Institute for Particle Astrophysics &; Cosmology, Depart- 
ment of Physics, and SLAC National Accelerator Laboratory Stan- 
ford University, Stanford, CA 94305 

3 TABASCO Fellow, Physics Department, University of Cali- 
fornia at Santa Barbara, 2233B Broida Hall, Santa Barbara, CA 
93106 

4 Fermi National Accelerator Laboratory, P.O. Box500, Batavia, 
IL 60510 

J Department of Physics, The University of Chicago, Chicago, 
IL 60637 

6 Kavli Institute for Cosmological Physics, The University of 
Chicago, Chicago, IL 60637 

7 Physics Department, University of Michigan, Ann Arbor, MI 
48109 

8 Astronomy Department, University of Michigan, Ann Arbor, 
MI 48109 

9 Michigan Center for Theoretical Physics, Ann Arbor, MI 48109 

10 Department of Astronomy and Astrophysics, The University 
of Chicago, Chicago, IL 60637 

11 University of California Observatories &; Department of As- 
tronomy, University of California, Santa Cruz, CA 95064 

12 Department of Physics & Astronomy, Northwestern Univer- 
sity, Evanston, IL 60208 

13 Brookhaven National Laboratory, Upton, NY 11973 

14 Department of Astronomy, The Ohio State University, Colum- 
bus, OH 43210 



et al. 2006; Tinker et al. 2008). Despite these successes, 
realizing the promise of cluster cosmology has proven dif- 
ficult. Indeed, a review of observational results from the 
past several years yields a plethora of studies where typ- 
ical uncertainties are estimated at the Acg ~ 0.05 — 0.10 
level despite a spread in central values that range from 
a s ps 0.65 to cr 8 ps 1.0 (Viana & Liddle 1996, 1999; Henry 
& Arnaud 1991; Henry 2000; Pierpaoli et al. 2001; Bor- 
gani et al. 2001; Seljak 2002; Viana et al. 2002; Schuecker 
et al. 2003; Allen et al. 2003; Bahcall et al. 2003; Bahcall 
& Bode 2003; Henry 2004; Voevodkin & Vikhlinin 2004; 
Rozo et al. 2007b; Gladders et al. 2007; Rines et al. 2007). 

The discrepancies among the various studies men- 
tioned above is a manifestation of the fundamental prob- 
lem confronting cluster abundance studies: theoretical 
predictions tell us how to compute the abundance of ha- 
los as a function of mass, but halo masses are not ob- 
servable. Consequently, we are forced to rely on observ- 
able quantities such as X-ray temperature, weak lensing 
shear, or other such signals, to estimate cluster masses. 
This reliance on observable mass tracers introduces sig- 
nificant systematic uncertainties in the analysis; indeed, 
this is typically the dominant source of error (e.g. Henry 
et al. 2008). 

There are two primary ways in which these difficul- 
ties can be addressed. One possibility is to reduce these 
systematic uncertainties through detailed follow-up ob- 
servations of relatively few clusters, an approach exem- 
plified in the work of Vikhlinin et al. (2008a). The second 
possibility is to use large cluster samples complemented 
with statistical properties of the clusters that are sensi- 
tive to mass to simultaneously fit for cosmology and the 
observable-mass relation of the cluster sample in ques- 
tion. Indeed, this is the basic idea behind the so called 
self-calibration approach, in which one uses the cluster- 
ing of clusters (Schuecker et al. 2003; Estrada et al. 2008) 
and cluster abundance data to derive cosmological con- 
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straints with no a-priori knowledge of the observable- 
mass relation (Hu 2003; Majumdar & Mohr 2004; Lima 

6 Hu 2004, 2005). There are, however, many other sta- 
tistical observablcs that correlate well with mass, such 
as the cluster-shear correlation function (Sheldon et al. 
2007), or even counts binned in multiple mass tracers 
(Cunha 2008). By including such data we can break 
the degeneracy between cosmology and the observable- 
mass relation, thereby obtaining tight cosmological con- 
straints while simultaneously fitting the observable-mass 
relation. 

In this work, we derive cosmological constraints from 
the SDSS maxBCG cluster sample (Koester et al. 2007a) 
and the statistical weak lensing mass measurement from 
Johnston et al. (2007b). We then compare our result to 
three state-of-the-art cluster abundance studies of X-ray 
selected cluster samples (Mantz et al. 2008; Henry et al. 
2008; Vikhlinin et al. 2008b) and demonstrate that our 
results are both consistent and competitive with these 
studies. This is the first time an optically selected catalog 
with masses estimated in a statistical way has produced 
constraints that are of comparable accuracy to the more 
traditional approach. 

The paper is organized as follows. Section 2 presents 
the data used in our study. Section 3 describes our anal- 
ysis, including the likelihood model and priors adopted 
in this work, and the way in which the analysis was im- 
plemented. Section 4 presents our main results, while 
sections 5.1, 5.2, and 5.3 discuss various sources of sys- 
tematic uncertainties. Section 6.1 compares our results 
to the most recent results from X-ray selected cluster 
samples. Section 6.2 investigates the implications of 
our results for dark energy, and Section 6.3 discusses 
the prospects for improving our cosmological constraints 
from the maxBCG cluster sample in the future. Section 

7 summarizes our main results and conclusions. Unless 
otherwise stated, all masses in this work are defined us- 
ing an overdensity A = 200 relative to the mean matter 
density of the universe. 

2. DATA 
2.1. MaxBCG Cluster Counts 

The maxBCG cluster catalog (Koester et al. 2007a) 
is an optically selected catalog drawn from 7,398 deg 2 
of DR4+ imaging data of the Sloan Digital Sky Sur- 
vey (SDSS). 15 The maxBCG algorithm exploits the 
tight E/S0 ridgeline of galaxies in color-magnitude space 
to identify spatial overdensities of bright red galaxies. 
The tightness of the color distribution of cluster galax- 
ies greatly suppresses the projection effects that have 
plagued optically selected cluster catalogs, and also al- 
lows for accurate photometric redshift estimates of the 
clusters (Az w 0.01). MaxBCG clusters are selected such 
that their photometric redshift estimates are in the range 
Zphoto £ [0.1,0.3], resulting in a nearly volume-limited 
catalog. A detailed discussion of the maxBCG cluster 
finding algorithm can be found in Koester et al. (2007b). 

We bin the maxBCG cluster sample in nine richness 
bins spanning the range -/V200 € [11, 120], corresponding 
roughly to M G [7 x lO 13 /^ 1 Mq, 1.2 x lO 15 /^ 1 M ]. 
Our richness measure N 2 oo is defined as the number of 

15 Wo write DR4+ as the catalog used a few hundred degrees of 
imaging beyond those released with DR4. 



TABLE 1 

Abundance of maxBCG 
Clusters 



Richness 


No. of Clusters 


11-14 


5167 


14-18 


2387 


19-23 


1504 


24-29 


765 


30-38 


533 


39-48 


230 


49-61 


134 


62-78 


59 


79-120 


31 
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Fig. 1. — Observed (diamonds) and modeled (solid line) cluster 
counts as a function of richness in our best-fit model described in 
Section 4. The model counts are computed using the best fit model 
detailed in Section 4, and are a good fit to the data. 

red-sequence galaxies within a scaled radius such that 
the average galaxy overdensity interior to that radius is 
200 times the mean galaxy density of the universe (see 
Koester et al. 2007a, for further details). The richness 
bins, and the number of clusters in each bin, are pre- 
sented in Table 1. There are an additional five clus- 
ters with richness N200 > 120. These five clusters have 
7V 2 oo = 126, 139, 156, 164, and 188, and are properly 
included in the analysis on an individual basis (see §3.1 
for details). 

Figure 1 shows the cluster counts corresponding to Ta- 
ble 1. Error bars between the various points are corre- 
lated. Also shown are the modeled counts from our best- 
fit model, detailed in Section 4. We show these model 
counts here for comparison purposes. 

2.2. MaxBCG Weak Lensing Masses 

Estimates of the mean mass of the maxBCG clusters 
as a function of richness are obtained through the weak 
lensing analysis described by Sheldon et al. (2007) and 
Johnston et al. (2007b). Briefly, Sheldon et al. (2007) 
binned the maxBCG cluster sample in richness bins as 
summarized in Table 2. Given a cluster in a specified 
richness bin, they use all cluster-galaxy pairs with the 
selected cluster as a lens to estimate the density contrast 
profile AS of the cluster. While these individual clus- 
ter profiles are very low signal-to-noise, averaging over 
all clusters within a richness bin allows one to obtain 
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TABLE 2 

Mean Mass of maxBCG Clusters 



Richness 


No. of Clusters 


(M 200!) ) [10 14 M Q ] 


12-17 


5651 


1.298 


18-25 


2269 


1.983 


26-40 


1021 


3.846 


41-70 


353 


5.475 


71+ 


55 


13.03 



Note. — Masses listed here are based on those 
quoted in Johnston et al. (2007b), rescalcd by the 
expected photometric redshift bias described in 
the text, and extrapolated to a matter overden- 
sity A = 200 from the A = 180 value quoted in 
Johnston et al. (2007b). The masses have also 
been rescaled to the cosmology that maximizes 
our likelihood function, (erg = 0.80, Q m = 0.28). 
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accurate estimates for the mean density contrast profile 
of maxBCG clusters function of richness. The re- 
sulting profiles arc fit using a halo model formalism to 
derive mean cluster masses by Johnston et al. (2007b). 
We then correct these masses upward by a factor of 1.18 
due to the expected photometric redshift bias due to the 
dilution of the lensing signal from galaxies that are in 
front of the cluster lenses, but whose photometric red- 
shift probability distribution extends past the cluster lens 
(see Mandelbaum et al. 2008b, for details). A very simi- 
lar but independent analysis has also been carried out by 
Mandelbaum et al. (2008a), and we use the comparison 
between the two independent analysis to set the system- 
atic error uncertainty of the weak lensing mass estimates 
(Rozo et al. 2008a). The final results of the weak lensing 
analysis summarized above are presented here in Table 
2. 16 Figure 2 shows the mean weak lensing masses from 
Table 2. Also shown are the mean masses computed us- 
ing the best-fit model detailed in Section 4. The richness 
binning of the weak lensing mass estimates differs from 
that of the abundance data because of the larger number 
of clusters necessary within each richness bin to obtain 
high S/N weak lensing measurements. 

3. ANALYSIS 

We employ a Bayesian approach for deriving cosmolog- 
ical constraints from the maxBCG cluster sample. We 
use only minimal priors placed on the parameters gov- 
erning the richness-mass relation, relying instead on the 
cluster abundance and weak lensing data to simultane- 
ously constrain cosmology and the richness-mass relation 
of the clusters. Details of the model, parameter priors, 
and implementation can be found below. 

3.1. Likelihood Model 

The observable vector x for our experiment is com- 
prised of: 

1. N% through Ag: the number of clusters in each of 
the nine richness bins defined in Table 1. 

2. (NM)i through (AM)s, the total mass contained 
in clusters in each of the five richness bins defined 

16 The number of clusters in Table 2 is larger than that reported 
in Johnston et al. (2007b) due to masking in the weak lensing mea- 
surements. This additional masking does not bias the recovered 
masses in any way. 



Fig. 2. — Mean weak lensing mass of maxBCG clusters as a 
function of richness. The diamonds with error bars correspond to 
our data, while the solid line shows the values predicted from our 
best-fit model (see Section 4 for details). We note the error bars 
are correlated, and the model is a good fit to the data. 

in Table 2, computed assuming Cl m = 0.27 and 
h = 0.71. 17 

We adopt a Gaussian likelihood model, which is fully 
specified by the mean and covariance matrix of our ob- 
servables. Expressions for these quantities as a function 
of model parameters are specified below. We also mul- 
tiply this Gaussian likelihood by a term that allows us 
to properly include the information contained in clusters 
with richness A200 > 120. In this richness range clus- 
ters are very rare and a Gaussian likelihood model is not 
justified. Instead we adopt a likelihood model where the 
probability of having a cluster of a particular richness 
A200 is binary (i.e. a Bernoulli distribution), with 

w^V™:!. a) 

Such a probability distribution is adequate so long as 
the probability of having two clusters of a given richness 
is infinitcsimally small. Note that given this binary prob- 
ability distribution, we have that the expectation value 
of the number of such clusters is simply (A(A2oo)) = P, 
and the likelihood is fully specified by the expectation 
value of our observable. We find that the likelihood of 
observing the particular richness distribution found for 
the maxBCG catalog for clusters of richness A200 > 120 
is 

C t ail= II 0- - ( N ( N 200)}) I] ( N ( N ™o)} ■ 
N{N 200 )=0 N(N 200 ) = 1 

(2) 

The first product is over all richness A 2 oo > 120 and 
no clusters in them, and the second product is over rich- 
ness bins which contain one cluster. The subscript tail 
reflects the fact that it is the likelihood of the tail of the 
abundance function. The final likelihood C = CoCtaa 
is the product of the Gaussian likelihood Lq described 
earlier and the likelihood of the abundance function tail. 

17 While Table 2 reports the masses after corrections, assum- 
ing Q m = 0.28, the actual input to our statistical analyis are the 
uncorrected masses from Johnston et al. (2007b), which assume 
n m = 0.27. 
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We note that the log-likelihood of the tail simplifies to 

A r 2oo>120 

]T {N(N 20O ))+hi{N(N 2O0 )).(3) 

N(N 200 ) = 1 

An identical result is obtained assuming only Poisson 
variations in the number of clusters for A^oo > 120. 

3.2. Expectation Values 

To fully specify our likelihood model we need to derive 
expressions for the mean and variance of our observables. 
The model adopted in this work is very similar in spirit 
to that of Rozo et al. (2007b), so we present here only a 
brief overview of the formalism. Interested readers can 
find a detailed discussion in Rozo et al. (2007b). 

We begin by considering the expected mean number 
of clusters in our sample. The number of halos within 
a rcdshift bin z 6 [z m i n , z max ] and within a mass range 
[M mm ,M m 

ax 

] is given by 



N = J dMdz ^M%^ M ^ 



(4) 



where dn/dM is the halo mass function, dV/dz is the co- 
moving volume per unit redshift, and ip(M) and <j>(z) are 
the mass and redshift binning functions, i.e. ip(M) = 1 if 
M is within the mass bin of interest and zero otherwise, 
and 4>{ z ) = 1 if z is within the redshift bin of interest, 
but is zero otherwise. 

In practice, we observe neither a cluster's mass nor its 
true redshift, but are forced to rely on the cluster rich- 
ness A^oo as a mass tracer and to employ a photometric 
rcdshift estimate. Let then P(N 2 qo\M) denote the prob- 
ability that a cluster of mass M has a richness A200 > and 
let P(z p hoto\z) denote the probability that a cluster at 
rcdshift z is assigned a photometric redshift z v h oto - The 
binning function ip is now a function of richness rather 
than mass so ^200) _= 1 for N 200 € 
Likewise, the redshift binning function is now a function 
of photometric redshift z p hoto- The total number of clus- 
ters in the maxBCG catalog becomes 



where 



dMdz^^-m-i) w), 



(V>|M>= / dN 200 P{N 200 \M)^{N 200 ) 



z) — dz p hoto P{Zphoto\z)4>(z p i lo to) 



(5) 



(6) 
(7) 



The quantity (ip\M) represents the probability that a 
halo of mass M falls within the richness bin defined by 
ip. We show these probabilities as a function of mass 
for each of the nine richness bins considered here in Fig- 
ure 3. To make the figure, we have set all relevant model 
parameters to their best-fit value detailed in Section 4. 

A similar argument allows us to write an expression 
for the expectation value for the total mass contained in 
clusters of a specified richness and redshift bin. This is 
given by 

(ATM) = I dMdz^^-M Q>\M) (4>\z) . (8) 




10" 10' 
Cluster Moss (M 200b , Solar Units) 
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Fig. 3. — Mass selection function of the maxBCG algorithm. The 
nine solid curves represent the probability that a halo of the corre- 
sponding mass falls within each of the nine richness bins described 
in Table 1. The dashed line is the sum of all the binning functions, 
and is the probability that a halo of a given mass is assigned a 
richness A^oo £ [H,120], i.e. it is the mass selection function of 
the maxBCG algorithm over this richness range. These binning 
functions are all estimated using our best-fit model parameters, 
which are detailed in Section 4. 

The notation NM reflects the fact that if M is the mean 
mass of the clusters of interest, the total mass contained 
in such clusters is NM where N is the total number of 
clusters in said bin. 

So far, our formulae adequately describe our exper- 
iment provided the weak lensing masses estimated by 
Johnston et al. (2007b) are fair estimates of the mean 
mass of the maxBCG clusters. In practice, there is an 
important systematic that needs to be properly incorpo- 
rated in our analysis, and which slightly modifies our ex- 
pression. We are referring to uncertainties in the photo- 
metric redshift estimates of the source galaxies employed 
in the weak lensing analysis. The main problem here 
is that the mean surface mass density profile S(i?) re- 
covered by the weak lensing analysis is proportional to 
1/ (S" 1 ), the average inverse critical surface density of 
all lens-source pairs employed in the analysis. We intro- 
duce an additional weak lensing bias parameter (3 such 
that if M true is the true mean mass of a set of clus- 
ters, the weak lensing mass estimate Mobs is given by 
M b s = P Mtrue- Consequently, our final expression for 
the mean weak lensing masses of the maxBCG clusters 



{NM) = (3 



dr> dV 

dMdZ jM^ M{m)m - (9) 



Priors on the parameter (3 are discussed in Section 3.4. 

3.3. Covariance Matrix 

There are multiple sources of statistical uncertainty in 
the data. These include: (1) Poisson fluctuations in the 
number of halos of a given mass, (2) variance in the mean 
overdensity of the survey volume, and (3) fluctuations in 
the number of clusters at fixed richness due to stochastic- 
ity of the richness mass relation. The covariance matrix 
of the observables is defined by the sum of the covariance 
matrices induced by each of the three sources of statisti- 
cal fluctuation just mentioned. A detailed derivation of 
the relevant formulae is presented in Rozo et al. (2007a). 
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Since this derivation generalizes trivially to include the 
mean mass as an additional observable — one needs only 
to introduce a mass weight in the formulae as appropriate 
— we will not repeat ourselves here. 

There is, however, one additional source of statisti- 
cal uncertainty that is not included in these calcula- 
tion, namely measurement error in the weak lensing 
masses. More specifically, uncertainties in the recovered 
weak lensing masses is dominated by shape noise in the 
source galaxies. This error was estimated by Sheldon 
et al. (2007) using jackknife resampling, and was prop- 
erly propagated into the computation of the weak lensing 
mass estimates by Johnston et al. (2007b). This error 
is added in quadrature to the diagonal elements of the 
covariancc matrix corresponding to the mean mass mea- 
surements. 

Finally, in addition to the errors summarized above, 
the covariance matrix is further modified due to system- 
atic uncertainties in the purity and completeness of the 
sample. The basic set up is this: if N trU e{N 2 oo) is the 
number of clusters one expects in the absence of systcm- 
atics, and N }, s (N 2 oo) is the actual observed number of 
clusters, one has 

N obs {N 200 ) = \(N 20 o)N true (N 200 ) 7 (10) 

where A is a factor close to unity that characterizes the 
purity and completeness systematics. If the sample is 
pure but incomplete, A is simply equal to the sample's 
completeness. For a complete but impure sample, A is 
one over the sample's purity. Note that, in general, A 
is itself a function of the cluster richness N 2 oo ■ In Bozo 
et al. (2007b), we estimated the purity and complete- 
ness of the maxBCG cluster sample at 95% or higher for 
A 2 oo > 11 (sec Figures 3 and 6 in that paper), suggesting 
A = 1.00 ± 0.05. Given that 

Var(A obs ) = Var(A)A 2 , ue + A 2 Var(A true ), (11) 

it follows that we can incorporate the impact of this nui- 
sance parameter by simply adding in quadrature the rela- 
tive uncertainty introduced by A to the covariance matrix 
estimated in the previous section. A similar argument 
holds for the total mass contained in clusters within each 
richness bin. That is, if M tr .„ e _is the true mean mass of 
clusters of richness A 20 o, and M b s is the observed mean 
mass, we expect 

(NM) obs = \{NM) true , (12) 

where A is a correction factor that accounts for the mass 
contribution of impurities in the sample. Unfortunately, 
it is impossible to know a priori what this factor A should 
be, even if we knew the correction factor A for cluster 
abundances. The reason is that false cluster detections 
will most certainly have a mass overdensity associated 
with them, just not that of a halo of the expected mass 
given the observed richness. Without a priori knowledge 
of this mass contribution, it is impossible to estimate 
the proper value of A. In the extreme case that all false 
detections have mass M, then the recovered value for 
NM will be biased by a factor A — A, which suggests 
adopting a fiducial value Var(A) = 0.05 2 to add to the 
diagonal matrix elements corresponding to the observed 
weak lensing masses. That is the approach we follow 
here. Throughout, we always set Var(A) = Var(A). 



3.4. Model Parameters and Priors 

Our analysis assumes a neutrino-less, flat ACDM cos- 
mology, and we fit for the values of <7s an d Q m - The 
Hubble parameter is held fixed at h = 0.7, and the tilt 
of the primordial power spectrum is set to n = 0.96 as 
per the latest WMAP results (Dunkley et al. 2008). The 
baryon density f2f,/i 2 is also held fixed at its WMAP5 
value fli,h 2 = 0.02273. Of these secondary parameters, 
the two that are most important are the Hubble constant 
and tilt of the primordial matter power spectrum (Rozo 
et al. 2004). Section 5.2.1 demonstrates our results are 
robust to marginalization over these additional parame- 
ters. 

The richness-mass relation P(A 2 oo|M) is assumed to 
be a log-normal of constant scatter. The mean log- 
richness at a given mass (\nN 2 Qo\M) is assumed to 
vary linearly with mass, resulting in two free parame- 
ters. We comment on possible deviations from linear- 
ity in Section 5.3.1. For the two parameters specify- 
ing the mean richness-mass relation we have chosen the 
value of (lnN 200 \M) at M = 1.3 x 10 14 M G and at 
Al = 1.3 x 10 15 Mq. These two are very nearly the val- 
ues of the mean mass for our lowest and highest richness 
bins, and therefore roughly bracket the range of masses 
probed in our analysis. The value of (In A 2 oo|M) at any 
other mass is computed through linear interpolation. We 
adopt flat priors on both of these parameters. 

The scatter in the richness-mass relation <Jn 200 \m is 
defined as the standard deviation of In A^oo at fixed M, 
tJ N 20 o\M = Var(ln A 200 |M). We assume that this quan- 
tity is a constant that does not scale with mass, and 
adopt a flat prior ct/v 2Q0 |m £ [0.1,1.5] for this parame- 
ter. We comment on possible deviations from constant 
scatter in Section 5.3.2. The minimum scatter allowed in 
our work (<Jn 200 \m = 0.1) corresponds to a 10% scatter, 
which is the predicted scatter for Yx in simulations. Yx 
is usually regarded as the X-ray mass tracer that is most 
tightly correlated with mass, so our prior on the scatter 
is simply the statement that richness estimates are less 
faithful mass tracers than Yx- 

We also place a prior on the converse scatter, that 
is, the scatter in mass at fixed richness a 2 r , Ar = 

Var(lnM|A 2 oo) at A 2 oo = 40. We emphasize that in our 
analysis the scatter om|at 2 oo ^ s considered an observable, 
not a parameter (the parameters is oat 200 |m)- The prob- 
ability distribution P(a^m 200 ) is taken directly from the 
analysis by Rozo et al. (2008a), and can be roughly sum- 
marized as om|tv 20 o = 0-45 ± 0.10 (1 a). This constraint 
is derived by demanding consistency between the ob- 
served Lx~N 2 qq relation of maxBCG clusters, the mass- 
richness relation of maxBCG clusters derived from weak 
lensing, and the Lx — M relation of clusters measured 
in the 400d survey (Vikhlinin et al. 2008b) . To compute 
the observed scatter o~m\n-, o as a function of our model 
parameters we directly compute the variance in log-mass 
for clusters in a richness bin A 2 oo € [38,42]. The vari- 
ance in In M due to the finite width of the bin is of order 
(1/40) 2 « 0.006, which is to be compared to the intrinsic 
variance « 0.45 2 « 0.2. Because the intrinsic variance is 
significantly larger than the variance due to using a fi- 
nite bin width, our results are not sensitive to the width 
of the bin used in the implementation of the prior. We 
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TABLE 3 



Parameter" 


Prior 6 


Importance 




[0.4,1.2] 


unrestrictive 




[0.05,0.95] 


unrcstrictive 


HniV200 Mi) 


flat 


unrestrictive 


(\nN 2 Qo\M 2 ) 


flat 


unrcstrictive 


a N 200 \M 


[0.1,1.5] 


unrestrictive 


P 


1.00 ±0.06; [0.5,1.5] 


restrictive 


a M\N 200 d 


Rozo et al. (2008a) 


restrictive 



a The masses Mi and M 2 are set to 1.3 X 10 14 M 
and 1.3 X 10 15 Mq respectively. 

b Priors of the form [a, b] mean the parameter in ques- 
tion is restricted to values within that range (flat 
prior). Priors of the form x = a ± 5a refer to a Gaus- 
sian prior of mean (x) = a and variance Var(x) = 
(5a) 2 . 

c Column specifies whether our results are sensitive 
to the assumed priors. We refer to a prior as restric- 
tive if our cosmological constraints are sensitive to the 
assumed prior, and unrestrictive otherwise. The only 
restrictive priors are that on the mass bias parame- 
ter /3 and the prior on the scatter in mass at fixed 
richness. 

d Note o"m|jv 2 oo * s n °t reau y a parameter in our anal- 
ysis but an observable that can be computed given the 
six parameters above. 

have explicitly checked that this is indeed the case. We 
have also checked that our results are insensitive to the 
location of the richness bin. That is, placing our prior 
on a M | Nwo at 7V 2 oo = 30 and iV 2 oo = 50 gives results 
that are nearly identical to those obtained with our fidu- 
cial -/V200 = 40 value. Finally, we note that in using the 
scatter measurement of Rozo et al. (2008a), who used an 
overdensity threshold of 500 relative to critical to define 
cluster masses, we are making the implicit assumption 
that the value of the current uncertainties in the scatter 
are much larger than any sensitivity to differences in the 
cluster mass definition. To address this concern, in Sec- 
tion 5.2.2 we discuss how the scatter prior impacts our 
results. 

The redshift selection function P(z p hoto\z) is assumed 
to be Gaussian with (z p hoto\z) = z and a(z p h oto \z) = 
0.008, as per the discussion in Koester et al. (2007a). 
We have explicitly checked that our results are not 
sensitive to our choice of parameters within the range 
8 {z P hoto\z) ~ 0.005 and 5a Zphoto \ z = 0.02, which encom- 
pass the uncertainties in the photometric redshift distri- 
bution of the maxBCG clusters (Koester et al. 2007a). 

Finally, we also adopt a prior on the weak lensing mass 
bias parameter, j3 = 1.0 ± 0.06, and allow it to vary over 
the range [0.5,1.5]. The width of our Gaussian prior is 
simply the mean difference between the Johnston et al. 
(2007b) masses (after correcting for photometric redshift 
bias) and those of Mandclbaum et al. (2008a) (for a more 
detailed discussion see Rozo et al. 2008a). 

The total number of parameters that are allowed to 
vary in our Monte Carlo Markov Chain (MCMC) is six: 
a s , dm, (lniV 2 oo|M) evaluated at M = 1.3 x 10 14 M 
and M = 1.3 x 10 15 Mq, aj^ %00 \M, and (3. We summarize 
the relevant priors in Table 3. 

3.5. Implementation 

We use the low baryon transfer functions of Eisenstein 
& Hu (1999) to estimate the linear matter power spec- 



trum. The halo mass function is computed using Tinker 
et al. (2008). We use a mass definition corresponding 
to a 200 overdensity with respect to the mean matter 
density of the universe, and adopt the Sheth-Tormen ex- 
pressions for the mass dependence of halo bias (Sheth & 
Tormen 2002) (this enters into our analysis only in the 
calculation of sample variance). The likelihood function 
is sampled using a Monte Carlo Markov Chain (MCMC) 
approach with a burn in of 22,000 points during which 
the covariancc matrix of the parameters is continually 
updated so as to provide an ideal sampling rate (Dunk- 
ley et al. 2005). We then run the chains for 10 5 points, 
and use the resulting outputs to estimate the 68% and 
95% likelihood contours in parameter space. For further 
details, we refer the reader to Rozo et al. (2007a). 

The one point that is worth discussing here is our cor- 
rections for the dependence of the recovered weak lens- 
ing masses on the assumptions about cosmology used for 
the measurements. Johnston et al. (2007b) quote halo 
masses at an overdensity of 180 relative to the mean 
background of the universe. Given that we use a density 
contrast of 200 relative to mean in order to compute the 
halo mass function, we must re-scalc the observed masses 
to our adopted mass definition. Moreover, the weak lens- 
ing analysis assumed Q m = 0.27. Given a different mat- 
ter density parameter Q m , the quoted mass will no longer 
correspond to an overdensity of 180, but to an overden- 
sity of 180(0. 27/Q m ). We explicitly apply this re-scaling 
to the observed weak lensing masses at each point in our 
MCMC. In practice, there is also an additional correc- 
tion due to the dependence of the lensing critical surface 
density E c on the matter density parameter O m , as well 
as small corrections due to systematic variations in halo 
concentration with mass. However, these corrections are 
expected to be small, and are fully degenerate with the 
mass bias parameter f3, so we do not include them here. 
The rescaling of the weak lensing masses is done using 
the fitting formulae in Hu & Kravtsov (2003). 

4. RESULTS 

Figure 4 presents the 68% and 95% confidence regions 
for each pair of parameters in our fiducial analysis de- 
scribed in §3. Plots along the diagonal show the proba- 
bility distributions of each quantity marginalized over the 
remaining parameters. Upper left plot showing the prob- 
ability distribution of the mass parameter (3 also shows 
the prior (3 = 1.00 ± 0.06 as a dashed curve. Our best 
fit model is summarized in Table 4, and is defined as 
the expectation value of all of our parameters. To test 
that our best fit model is a good model to the data, we 
performed 10 4 Monte Carlo realizations of our best fit 
model, and evaluated the likelihood function for each of 
these realizations. Setting (ln£) = 0, from our Monte 
Carlo realizations we find ln£ = 0.0 ± 6.9, which is to 
be compared to the data likelihood ln£ = —5.2. The 
data likelihood is therefore consistent with our model, 
demonstrating the model is statistically a good fit. 

In the discussion that follows, we restrict ourselves 
to the subset of plots which we find most interesting. 
Throughout, unless otherwise noted we summarize con- 
straints on a parameter p by writing p = p + a p where 
p and a p are the mean and standard deviation of the 
likelihood distribution for p marginalized over all other 
parameters. We use this convention even when the like- 
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show 68% and 95% confidence regions. Plots along the diagonal show the probability distributions for each quantity marginalized over the 
remaining parameters. The probability distribution for the mass bias parameter fi also shows the prior f3 = 1.00 ± 0.06 assumed in the 
analysis. 



TABLE 4 
Best Fit Model 



Parameter 


maxBCG 


maxBCG+WMAP5 b 


as 


0.804 ± 0.073 


0.807 ±0.020 




0.281 ±0.066 


0.269 ±0.018 


(lnTVaoolAii) 


2.47 ±0.10 


2.48 ±0.10 


(lniV200 M2} 


4.21 ±0.19 


4.21 ±0.13 


a N 200 \M 


0.357 ±0.073 


0.348 ± 0.071 




1.016 ±0.060 


1.013 ±0.059 



a The masses Mi and M 2 are set to 1.3 X 10 14 M 
and 1.3 X 10 15 Mq respectively. 

b Thcse values are obtained by including the WAMP5 
prior <T 8 (fW0.25)-°' 312 = 0.790 ± 0.024. See Section 
4.3 for details. 



lihood function is obviously not Gaussian. 

4.1. Cosmological Constraints and Comparison to 
WMAP 

The solid curves in Figure 5 show the 68% and 95% 
confidence regions from our analysis. The "thin" axis 
of our error ellipse corresponds to os(n m /0.25) OA1 = 



0.832 ± 0.033. 18 The constraints on each of the in- 
dividual parameters are as = 0.80 ± 0.07 and fi m = 
0.28 ± 0.07. The marginalized likelihood can be rea- 
sonably approximated by a log-normal distribution with 
lnO TO = -1.313 ± 0.183, (lna 8 ) = -0.219 ± 0.081, 
and a correlation coefficient between lnf2 m and lnag 
r = —0.899. Also shown in Figure 5 as dashed curves 
are the corresponding regions from the WMAP 5-year 
results (Dunklcy et al. 2008). Our results are consistent 
with WMAP5. Combining the two experiments results 
in the inner filled ellipses, given by as = 0.807±0.020 and 
f2 m = 0.265 ± 0.016, with nearly no covariance between 
the two parameters (r = 0.008). These joint constraints 
on as and TO represent nearly a factor of two improve- 
ment relative to the constraints from WMAP alone. 

The shape of the confidence region is easy to interpret: 
since the number of massive clusters increases with both 
as and f2 m , in order to hold the cluster abundance fixed 
at its observed value any increase in as must be compen- 
sated by a decrease in O m , implying that a product of the 
form as^Zi must be held fixed. The specific value of 7 de- 

18 The exponent 0.41 is obtained by estimating the covariance 
matrix of In as and In Q m , and finding the best constrained cigen- 
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Fig. 5. — Constraints on the ag — Q m plane from maxBCG and 
WMAP5 for a flat ACDM cosmology. Contours show the 68% and 
95% confidence regions for maxBCG (solid), WMAP5 (dashed), 
and the combined results (filled ellipses). The thin axis of the 
maxBCG-only ellipse corresponds to cr 8 (n m /0.25) 41 = 0.832 ± 
0.033. The joint constraints are as = 0.807 ± 0.020 and n m = 
0.265 ± 0.016 (one-sigma errors). 
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Halo Mass M, 



Fig. 6. — Halo mass function for two different cosmologies satis- 
fying the maxBCG constraint o- 8 (^m/0.25) - 41 = 0.832. The mass 
functions arc weighted by the volume probed by the maxBCG cat- 
alog (computed assuming Q m = 0.265), and by the mass selection 
function shown in Figure 3. The maxBCG normalization condi- 
tion ijg(f2 m /0.25) ' 41 = 0.832 results in a fixed halo abundance at 
a mass scale M = 3.6 X 10 14 Mq. The dotted line at the top marks 
the mass scale at which the mean of the richness-mass relation is 
best constrained in our fiducial analysis. 



pends on the mass scale that is best constrained from the 
data. The particular degeneracy recovered by our anal- 
ysis corresponds to a mass scale M = 3.6 x 10 14 Mq, 
which is about what we would expect (i.e. roughly half 
way between the lowest and highest masses probed by our 
data). Figure 6 illustrates this argument by showing the 
Tinker et al. (2008) halo mass function weighted by the 
mass selection function from Figure 3 for two different 
cosmologies: a low eg (high f2 m ) cosmology, and a high 
cr 8 (low Q m ) cosmology, where the product er 8 ^m 41 has 
been held fixed to our best-fit value. We will refer back 
to Figure 6 multiple times in the following discussion. 

4.2. Constraints on the Richness-Mass Relation 



In our analysis, we parameterized the richness-mass re- 
lation in terms of its scatter, and the value of the mean 
(In N 200 \M} at two mass scales, M 1 = 1.3 x 10 14 M Q and 
M2 = 1.3 x 10 15 Mq. We now re-parameterize this rela- 
tion in terms of an amplitude and slope for (lniV2oo|M), 
selecting as the pivot point the mass scale at which the 
uncertainty in (\nN2Q0\M) is minimized. We write then 

(ln/V 20 o|M) = A + aQnM - lnM plvot ). (13) 

We find the error on the amplitude parameter is mini- 
mized for M p i VOt = 1.09 x 10 14 , which agrees well with the 
peak in the mass distribution of our clusters as shown in 
Figure 6. In what follows, we discuss only constraints on 
the richness mass relation assuming this parameteriza- 
tion. A discussion of possible curvature in the richness- 
mass relation and/or mass scaling of its scatter is rele- 
gated to §5.3. 

Figure 7 summarizes our constraints on the richness- 
mass relation after marginalizing over all other parame- 
ters. The best-fit values for each of the parameters are 
A = 2.34 ± 0.10, a = 0.757 ± 0.066, and o-jv 200 |m = 
0.357 ±0.073. Note that for a pure power-law abundance 
function, one expects ctjv 200 |j\/ = o<-o~m\n WO i in accordance 
with our result. 

Of these results, the constraints on the slope and scat- 
ter of the richness-mass relation are particularly worth 
noting. First, it is clear that the naive scaling -/V200 oc M 
is not satisfied, with the slope of the richness-mass rela- 
tion being significantly smaller than unity. Second, the 
recovered scatter (7jv 20Q |m = 0.357 ± 0.073 is larger than 
the Poisson value cr/v, 00 |M ~ 0.2 that one might naively 
expect for clusters with -/V200 ~ 30 galaxies, which is the 
typical richness of clusters at the mass scale where mass 
function is best constrained. 

Interpreting these results in terms of standard halo oc- 
cupation model parameters requires care. The maxBCG 
richness is known to suffer from various sources of sys- 
tematics including miscentering of clusters (Johnston 
et al. 2007a) and color off-sets in the richness estimates 
(Rozo et al. 2008b), both of which will impact the recov- 
ered richness-mass relation at some level. Moreover, any 
richness estimate will suffer to some extent from projec- 
tion effects (Cohn et al. 2007), and discrepancies between 
assigned cluster radii and the standard mass-overdensity 
definitions used for halos. Disentangling the various con- 
tributions of each of these different sources of scatter to 
the total variance of the richness-mass relation is beyond 
the scope of this paper, and will not be considered further 
here. 

Figure 7 also shows that the amplitude of the richness- 
mass relation is anti-correlated with the scatter. This is 
not surprising: at fixed cluster abundance, and given 
a fixed mass function, models with a high amplitude 
of the richness-mass relation result in halos that tend 
to be very rich. This means that the number of lower 
mass halos that scatter into higher richness must be 
low, or otherwise the abundance of clusters will be over- 
predicted. Consequently, high amplitude models must 
have low scatter, leading to an anti-correlation between 
the two parameters. 

4.3. Degeneracies Between Cosmology and the 
Richness-Mass Relation 
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Fig. 7. — Constraints on parameters of the richness-mass re- 
lation. Countours indicate the 68% and 95% confidence regions; 
diagonal histograms indicate the probability distribution for each 
parameter, marginalized over the remaining parameters. The am- 
plitude and slope parameters define the mean of the richness-mass 
relation as per Eq. 13. The pivot point of the relation occurs at 
M = 1.15 x 10 14 M e . 

Figure 4 shows that the most significant correlation be- 
tween cosmology and our fiducial richness-mass relation 
parameters is that between as and (In A200I-M2) where 
M 2 is our higher reference mass M 2 = 1.3 x 10 15 M Q . 
Because the pivot point for the mean of the richness- 
mass relation is so close to our original low mass reference 
scale Ah = 1.3 x 10 14 M Q used to define (In jV 20 o|-M"), it 
follows that M2 must be closely related to a, the slope 
of the richness-mass relation. We thus expect a strong 
degeneracy between as and a (see also Rozo et al. 2004). 

Figure 8 shows that this is indeed the case. We can 
understand the origin of this anti-correlation by investi- 
gating Figure 6. We have seen that the data fixes the 
amplitude of the halo abundance at M — 3.9 x 10 14 Mq. 
At the high mass end, however, the expected abundance 
of massive halos varies rapidly with a&. Low as models 
result in fewer massive halos, so high richness clusters 
will have relatively lower masses. That is, richness must 
increase steeply with mass, and hence a must be high, 
explaining the anti-correlation between as and a. 

Figure 8 also demonstrates how these constraints are 
improved when we include a WMAP five-year data prior 
cr 8 (fi m /0.25)-°- 312 = 0.790 ± 0.024. This prior corre- 
sponds to the error along the thin direction of the WMAP 
error ellipse. Since WMAP data breaks the as — Q m de- 
generacy in the data, including the WMAP prior pro- 
duces a tight constraint in the as — a plane. The new 
marginalized uncertainty in the slope of the richness- 
mass relation is a = 0.752 ± 0.024, significantly smaller 
than unity. 

5. SYSTEMATIC ERRORS 

We now consider the impact of three varieties of sys- 
tematic errors on our analysis. Section 5.1 investigates 
observational systematics, Section 5.2 investigates sys- 
tcmatics due to our assumed priors, and Section 5.3 in- 
vestigates systematics due to the parameterization of the 
richness-mass relation. 

5.1. Observational Systematics 

In this section, we study how observational systemat- 
ics affect the recovered cosmological constraints from our 
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Fig. 8. — Confidence regions in the erg — a plane. Solid ellipses 
show the 68% and 95% likelihood regions. The tight correlation 
between erg and a, the slope of the richness-mass relatio, can be 
understood on the basis of Figure 6: a low as implies few massive 
halos, so to avoid under-predicting the abundance of rich clusters, 
galaxies must preferentially live in lower mass halos, resulting in 
a more rapidly rising richness-mass relation (i.e. higher slope). 
This degeneracy is broken upon inclusion of the WMAP five-year 
constraint rj 8 (Q m /0.25)-°- 312 = 0.790 ± 0.024 as an additional 
prior, as illustrated by the inner dashed ellipses in the Figure. The 
corresponding constraint on the slope of the richness-mass relation 
is a = 0.752 ±0.024. 

analysis. We consider two such systematics: one, the im- 
pact of purity and completeness, and two, the impact of 
possible biases in the weak lensing mass estimates of the 
maxBCG clusters. We do not discuss uncertainties in the 
photometric redshifts for clusters at any length since, as 
discussed in Section 3.4, they are found to be negligi- 
ble. This is not surprising, as the maxBCG photometric 
rcdshift estimates are extremely accurate (a z « 0.008, 
Koester et al. 2007a). 

5.1.1. The Impact of Purity and Completeness 

Figure 9 compares the cosmological constraints ob- 
tained assuming perfect purity and completeness with 
those obtained assuming a 5% uncertainty in these quan- 
tities. While non-negligible, the 5% uncertainty in the 
completeness and purity function of the maxBCG cata- 
log is far from the dominant source of uncertainty in our 
analysis. Moreover, this uncertainty elongates the error 
ellipse along its unconstrained direction, but has a min- 
imal impact on the best constrained combination of as 
and £l m : A(78(^in/0.25) 41 = 0.033 in our fiducial anal- 
ysis, while Ac78(^m/0.25) - 41 = 0.029 assuming perfect 
purity and completeness, a mere 10% difference. 

It is easy to understand why a 5% uncertainty in the 
purity and completeness has a minimal impact in our re- 
sults. For A200 ^ 25, the statistical uncertainties in the 
cluster abundances are larger than the 5% uncertainty in 
the counts from purity and completeness. Since the best 
constrained combination of cosmological parameters is 
driven primarily by high mass clusters, a 5% uncertainty 
in the purity and completeness functions has little im- 
pact on this parameter combination. How far the error 
ellipse extends along the degeneracy, however, is primar- 
ily driven by the observational constraints on the low end 
of the halo mass function (see Figure 6). Consequently, 
the 5% systematic uncertainty in the low richness cluster 
counts elongates the error ellipse along its major axis. 

We conclude that for the expected level of purity and 
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Fig. 9. — Effect of purity, completeness, and mass bias on param- 
eter constraints. Plot shows 68% confidence regions assuming per- 
fect purity and completeness (dashed), increasing the width of the 
weak lensing mass bias prior from f3 = 1.00±0.05 to j3 = 1.00±0.12 
(dotted), and for our fiducial analysis (solid). We find the uncer- 
tainty in purity and completeness has a minimal impact on the best 
constrained combination of the erg and Q m parameter, and there- 
fore on the constraints from a joint maxBCG + WMAP5 analysis. 
The same is not true of the weak lensing mass bias parameter. Un- 
certainties in the maxBCG cluster masses are the dominant source 
of systematic in our current analysis, and increase the uncertainty 
of the parameter combination ag (r2 m /0.25) 0,4 by 45% (see Figures 
10 and 11 for comparison). 

completeness of the maxBCG cluster sample, our cosmo- 
logical constraints are robust to these systematics. 

5.1.2. Systematic Uncertainties of the Weak Lensing 
Mass Estimates 

In Section 2.2, we discussed that the weak lensing 
masses of Johnston et al. (2007b) were boosted by a fac- 
tor of 1.18 to account for biases arising from scatter in 
the photometric redshift estimates (Mandclbaum et al. 
2008b). Even with such a boost, the Johnston et al. 
(2007b) and the Mandelbaum et al. (2008a) mass esti- 
mates were not consistent, which led us in Section 3.2 to 
introduce a mass bias parameter /3 that uniformly scales 
all masses by the same amount in order to account for 
any remaining biases. We now wish to explore how ro- 
bust our results are to our estimate of this systematic 
uncertainty. 

Figure 9 illustrates what happens if we repeat our fidu- 
cial analysis while doubling the width of the prior of p 
from /3 from p = 1.00 ± 0.06 to /? = 1.00 ± 0.12. We 
find that the wider P prior significantly increases the 
uncertainty in the parameter combination o~g(£l m /0.25) 
from Aer 8 (ft m /0.25) a41 = 0.033 to Acr 8 (f2 m /0.25) - 41 = 
0.045, corresponding to a 36% increase of the error 
bar. Using this new, wider prior, we find that the 
joint maxBCG + WMAP 5-year likelihood result in 
the cosmological constraints as = 0.802 ± 0.023 and 
O ro = 0.261 ± 0.019, which constitute a ~ 15% increase 
in the uncertainty of each of these parameters respec- 
tively. Even with this wider prior, however, adding the 
maxBCG constraint to the WMAP5 result improves the 
final cosmological constraints on as and f2 TO by a fac- 
tor of 1.6 relative to those obtained using WMAP data 
alone. 

We can understand the impact of the mass bias pa- 
rameter on our cosmological constraints using Figure 6. 



A wider prior on P implies that the mass scale of the 
maxBCG clusters is more uncertain, so the mass at which 
the cluster abundance is best constrained, i.e. the point 
at which the two curves in Figure 6 cross each other, is 
more uncertain. Consequently, the cluster normalization 
constraint agf^ 41 is weakened. The error along the long 
direction of the error ellipse does not change because the 
width of the mass range probed by the maxBCG clusters 
is largely independent of an overall mass bias. 

One of the curious results that we have found in our 
study of the mass bias parameter P is that the prior and 
posterior distributions of this parameter are different. In 
particular, we find that given the priors P = 1.00 ± 0.06 
and P = 1.00 ± 0.12, the posterior distributions for P are 
P = 1.02 ±0.06 and = 1.06 ±0.12 respectively. Indeed, 
this explains why the error ellipse for our wider prior 
is displaced to the left of that of our fiducial analysis: 
the shift in P corresponds to a change in the mass scale, 
which has to be compensated by a change in the matter 
density parameter fl m . 

We conclude that the uncertainty in the weak lensing 
mass estimates of the maxBCG clusters is an important 
source of systematic uncertainty in our analysis. In fact, 
it is the dominant source of systematic uncertainty in 
our analysis. We have explicitly considered the impact 
of photometric redshift estimates for source galaxies as 
the source of this uncertainty, but other biases to the 
lensing masses — for example if the fraction of miscen- 
tered clusters was over- or under-estimated by Johnston 
et al. (2007b) — would affect our results in a similar way. 

5.2. Prior-Driven Systematics 

Our analysis makes use of two important priors: that 
the only two cosmological parameters of interest are as 
and Q m , and that the scatter in the richness-mass rela- 
tion can be determined from X-ray studies as discussed 
in Rozo et al. (2008a). Here, we discuss how our results 
change if these priors are relaxed. 

5.2.1. Cosmological Priors 

After as and f2 m , cluster abundance studies arc most 
sensitive to the Hubble parameter h and the tilt n of 
the primordial power spectrum. In Figure 10, we illus- 
trate how the constraints on the as — Q m plane are af- 
fected upon marginalization over h and n using Gaus- 
sian priors h = 0.7 ± 0.1 and n = 0.96 ± 0.05. As we 
can see, marginalizing over the Hubble parameter and 
the tilt of the power spectrum elongates the error el- 
lipse, but it does not make it wider. Thus, the combi- 
nation og^m 41 remains tightly constrained, and a joint 
maxBCG and WMAP 5-year data analysis is robust to 
the details of the priors used for h and n when estimating 
the maxBCG likelihood function. We also investigated 
whether a non-zero neutrino mass could significantly af- 
fect our results. Using a prior m u < 1 eV, we find that 
massive neutrinos do not significantly affect our constrain 
on er 8 (f2 m /0.25) - 41 . We conclude that holding the Hub- 
ble parameter and the tilt of the power spectrum fixed 
does not result in systematic uncertainties in the joint 
maxBCG + WMAP 5-year data analysis. 

5.2.2. The Impact of the Scatter Prior 

In Rozo et al. (2008a), we derived an empirical con- 
straint on the scatter of the richness-mass relation by 
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Fig. 10. — Effect of relaxing additional cosmological parame- 
ter priors on erg and Q m constraints. Lines show 68% confidence 
regions for our fiducial analysis (solid), after marginalizing over 
Gaussian priors h = 0.7 ± 0.1 and n = 0.96 ± 0.05 (dashed), and 
using only a fiat prior on the scatter in mass at fixed richness 
a M\N 200 e [0.1,1.5] (dotted). We find that holding the Hubble 
parameter and power spectrum index fixed does not bias nor arti- 
ficially tighten our constraint on crgf!^ 41 . The scatter prior from 
Rozo et al. (2008a) on &m\n 2 oo em pl°y e d m our fiducial analysis is 
found to have a significant impact on our data. More specifically, 
dropping this scatter prior increases the error along the short axis 
of our error ellipses by 36%. We have also explored whether mas- 
sive neutrinos significantly impact our constraint on crgC^j 41 , and 
find that for neutrino masses m v < 1 eV there is no degradation 
of the error. 



demanding consistency between X-ray, weak lensing, and 
cluster abundance data. The recovered scatter, however, 
characterized the richness-mass relation using a mass 
that was defined using an overdensity of 500 relative to 
the critical density of the universe. In this analysis, we 
use a density threshold of 200 relative to mean, so the 
use of the X-ray derived scatter prior is justified only 
if the scatter in the mass scaling between the two over- 
density thresholds is not the dominant source of scatter. 
While we fully expect this assumption to hold, we have 
repeated our analysis without use of the scatter prior in 
order to cross-check our results. 

Figure 10 summarizes our results. We find that our 
scatter prior tightens the error ellipse along both its short 
and long axis. This is as expected: without the scatter 
prior, the mass scale of the maxBCG clusters becomes 
less constrained, and consequently the halo mass func- 
tion is less tightly constrained at all scales. The best 
constrained combination of o% and Q m when dropping 
the Rozo et al. (2008a) prior on the scatter in the mass- 
richness relation is cx 8 (£l m /0.25) - 48 = 0.841±0.045. This 
value represents a 36% increase in uncertainty relative to 
our fiducial analysis. The joint maxBCG + WMAP5 
constraints in this case are as = 0.805 ± 0.021 and 
n m = 0.264 ±0.017. 

Not surprisingly, prior knowledge of the scatter of the 
mass-richness relation can significantly enhance the con- 
straining power of the maxBCG data set. Neverthe- 
less, even without prior knowledge in the scatter the 
joint maxBCG+WMAP constraints improve upon the 
WMAP values by a factor of 1.7. 

5.3. Parameterization Systematica 



One of the most important systcmatics that need to be 
addressed in studies where the observable-mass relation 
is parameterized in some simple way is how to assess the 
robustness of the results to changes in the parameteri- 
zation of the observable-mass relation. Here, we have 
assumed that the richness-mass relation P(N2oo\M) i s 
a log-normal of constant scatter and that (lniV2oo|M} 
varies linearly with InM. We now investigate how our 
results change if we relax some of these assumptions. 

5.3.1. Curvature in the Mean Richness-Mass Relation 

To investigate the impact of curvature in the mass 
richness relation, we assume (lniV2oo|M) is a piecewise 
linear function. We first specify (lniV2oo|M) at three 
mass scales Mi, M 2 , and M3, and define the value of 
(ln/V2oo|M) at every other mass through linear inter- 
polation in log-space. We set the minimum and max- 
imum reference masses to the same values as before, 
Mi = 1.3 x 10 14 M 0! and M 3 = 1.3 x 10 14 M ? . The 
intermediate reference mass is set to the geometric aver- 
age of these two masses, lnM2 = 0.5(ln Mi + I11M3), or 
M2 = 3.66 x 10 14 Mq. Note this mass scale is very nearly 
the same as the mass at which the halo mass function is 
best constrained. 

Figure 11 shows how our cosmological constraints 
change with the introduction of mass dependence on the 
slope of the mean richness-mass relation (I11./V200IM). 
We find that the thin axis of the error ellipse is not 
significantly affected by this more flexible parameteriza- 
tion, while the long axis of the error ellipse is somewhat 
lengthened. This is as expected: the high mass end of 
the halo mass function is only sensitive to how richness 
varies with mass for large M, and in this regime the more 
flexible parameterization does not introduce significantly 
more freedom. Thus, our data will tightly constrain the 
high mass end of the halo mass function just as well 
as did before, leading to no degradation in the error of 
ofcfi^ 41 ■ Once the high mass end of the richness-mass re- 
lation has been fixed, however, introducing curvature in 
(In /V200IM) dilutes the information contained in the low 
mass end of the halo mass function, thereby increasing 
the error ellipse along its long axis. Note the robustness 
of the <78^m 41 constraint also implies that the constraints 
of a joint maxBCG + WMAP 5 analysis are not signifi- 
cantly affected by our choice of parameterization. 

Irrespective of the impact our new parameterization 
of (ln7V2oo|M) has on our cosmological constraints, it is 
fair to ask whether or not there is significant evidence 
for curvature of the mean richness-mass relation. Us- 
ing a maximum likelihood ratio test, we find that the 
increase in likelihood due to curvature in the richness- 
mass relation is significant at the 50% level, less than la. 
Thus, there is no evidence for curvature in the richness- 
mass relation. We have also explicitly confirmed that the 
slopes of the low and high mass end of the richness-mass 
relation are consistent with each other. Indeed, we find 



d 2 (In N 200 \M} 



d In M 2 



= 0.05 ±0.07 

Af=3.66xl0 14 M Q 

where we have assumed 
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Fig. 11. — Effect of relaxing assumptions about the richness- 
mass relation on ag-Q m constraints. Contours show 68% con- 
fidence limits for our fiducial analysis (solid curve), assuming 
(In A^oo \M) is a piece- wise linear function (dashed) , and allowing 
&N 200 \M t° vary linearly with mass (dotted). Giving additional 
freedom to the richness-mass relation has a minimal impact on 
our constraint on agO.^ 41 . Moreover, using a likelihood ratio test 
we find that there is no evidence in the data for curvature of the 
richness-mass relation, nor for a scatter that varies with mass. We 
conclude that our parameterization of (In N20o\M) and Vn 200 \m 
do not introduce any significant systematics in our analysis. 

5.3.2. Scaling of the Scatter in the Richness-Mass 
Relation with Mass 

We now investigate whether allowing the scatter of the 
richness-mass relation to vary with mass has a signifi- 
cant impact on our cosmological parameters. For these 
purposes, we allow the scatter (Tn wo \m to vary linearly 
with InM, and parameterize it by specifying its val- 
ues at the reference masses M\ = 1.3 x 10 14 Af Q and 
M 2 = 1.3 x 10 15 M Q . The value of a Nmo \ M at any other 
mass is obtained through linear interpolation. 

Figure 11 compares the cosmological constraints we 
obtain with our new model to those of our fiducial anal- 
ysis with constant scatter. Once again, we find that the 
"thin" axis of the error ellipse is not significantly affected 
by the new more flexible parameterization, while the long 
axis is slightly elongated. The interpretation of these re- 
sults is the same as those of §5.3.1. We have tested for 
evidence of scaling of the scatter in the richness-mass re- 
lation with halo mass using a likelihood ratio test. The 
increase in likelihood due to a linearly varying scatter is 
significant at the 39% level, implying there is no evidence 
of mass dependence in the scatter of the richness-mass 
relation in the data. We have also explicitly confirmed 
that the scatter at the low and high mass ends probed 
by the maxBCG cluster sample are consistent with each 
other. Indeed, our constraint on the slope of the mass 
dependence of the scatter in the richness-mass relation 
is 
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where we assumed 
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We note the velocity dispersion analysis in Becker et al. 
(see 2007) points towards some mass dependence in the 



Fig. 12. — Effect of the minimum richness on crg-f2 m constraints. 
Coutours show 68% confidence regions for our fiducial analysis 
(solid curve), an analysis where we include an additional richness 
bin, v = 9 — 10, at the low end of the richness function (dashed), 
and an analysis where we drop the lowest richness bins considered 
in our fiducial analysis. Our cosmological constraints arc consis- 
tent for all these analyses. We have also found that removing the 
most massive clusters from our analysis has minimal impact on our 
cosmological constraints. 

scatter of the mass-richness relation, though part of this 
discrepancy is likely due to miscentering systematics (see 
Rozo et al. 2008a, for details). We are now in the process 
of reanalyzing the velocity dispersion data updating both 
our treatment of systematics, and substantially increas- 
ing the sample of spcctroscopically sampled galaxies, so 
we defer a detailed discussion of these results to a future 
paper. 

We conclude that our parameterization of the mean 
and scatter of the richness-mass relation does not intro- 
duce systematic errors in our analysis. 

5.3.3. Richness Range Considered 

We have tested whether there is cosmological infor- 
mation in the richness range -/V200 > 120 by running 
MCMCs both with and without the contribution of these 
clusters to the likelihood function. We find that these two 
analyses yield nearly identical results. We have also ex- 
plicitly confirmed that our results are robust to the low- 
est richness bin employed in the analysis. As we might 
expect, removing the lowest richness bin increases our 
uncertainties along the long axis of the error ellipse as 
shown in Figure 12. We also investigate adding a new 
lowest richness bin, consisting of clusters in with A^oo = 
9-10, as well as the mean mass for clusters in the range 
^200 = 9-11. This analysis rotates the error ellipse very 
slightly compared to our fiducial analysis, but does not 
significantly affect our results. 

6. DISCUSSION 

6.1. Comparison to Other Work 

The main point of this section is to demonstrate two 
points: 

1. The cosmological constraints from the maxBCG 
cluster catalog are competitive with the state of 
the art constraints derived from low redshift X-ray 
selected cluster samples. 
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2. Despite the markedly different analyses and sources 
of systematic uncertainty, the cluster abundance 
constraints from the maxBCG cluster sample are 
in excellent agreement with those of X-ray selected 
samples. This demonstrates the robustness of clus- 
ter abundance studies as a tool of precision cosmol- 
ogy. 

Given our goal, in this section we focus exclusively on 
the most recent cosmological constraints derived from 
low redshift X-ray cluster samples. In particular, we ex- 
plicitly consider only three works: Mantz et al. (2008), 
who worked with the X-ray luminosity function, Henry 
et al. (2008), who worked with the X-ray temperature 
function, and Vikhlinin et al. (2008b), who estimated 
the low redshift halo mass function using the 400d X-ray 
survey (Burenin et al. 2007) with mass estimates based 
on Yx (Kravtsov et al. 2006). These three papers are 
the most recent analyses of X-ray selected cluster sam- 
ples, and all recover tight cosmological constraints that 
are in excellent agreement with one another, while care- 
fully accounting for the relevant systematics for each of 
their analyses. 

Now, as we have discussed in previous sections, the 
main result from low redshift cluster abundance stud- 
ies is a tight constraint on the value of og^m where for 
maxBCG clusters 7 = 0.41. Other cluster samples, how- 
ever, will have slightly different values of 7, which brings 
up the question of how can we fairly compare these var- 
ious constraints. One way would be to simply quote the 
percent uncertainty in the relevant os^m combination. 
However, we would like to have a clear graphical repre- 
sentation of this result. We have chosen to do this by 
plotting the 68% confidence regions of a simplified ver- 
sion of a joint cluster abundance + WMAP5 analysis 
assuming a neutrino-less flat ACDM cosmology. We pro- 
ceed as follows: given a cluster abundance experiment, 
we consider only the constraint on a^^Zn disregarding all 
other cosmological information. We then add a WMAP 
5 prior cr 8 (f2 m /0.25)-°- 312 = 0.790 ± 0.024, which corre- 
sponds to the thin axis of the error WMAP5 error ellipse 
in the o~s — Q m plane, and we compute the corresponding 
68% confidence regions in the a% — f2 m plane. 

The result of this exercise is shown in Figure 13. The 
specific constraints from each of the works considered 
here are presented in Table 5. The agreement among 
the different analyses is excellent despite the tight error 
bars and the different sources of systematic uncertain- 
ties. This agreement clearly demonstrates not only that 
optically selected cluster samples can produce cosmolog- 
ical constraints that are competitive with those of X-ray 
selected cluster samples, but also that systematic uncer- 
tainties have been properly estimated. 

6.2. Low Redshift Cluster Abundances and The 
Equation of State of Dark Energy 

Detailed analyses exploring how cluster abundances 
help improve dark energy constraints have been pre- 
sented by previous groups, most recently by Mantz et al. 
(2008) and Vikhlinin et al. (2008b). Rather than dupli- 
cating their work, in this section we opt for performing a 
simple analysis that captures the essential physics behind 
the Mantz et al. (2008) and Vikhlinin et al. (2008b) re- 
sults, which helps illustrate exactly why and how clusters 
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Fig. 13. — Comparion of optical and X-ray cluster abundance 
constraints on erg — Q m . Contours show 68% confidence regions 
for a joint WMAP5 and cluster abundance analysis assuming a 
flat ACDM cosmology. In addition to our results (filled ellipse), 
we consider the latest cluster abundance constraints from the low 
redshift cluster luminosity function (dashed Mantz et al. 2008), 
temperature function (dash-dot Henry et al. 2008), and mass func- 
tion as estimated with Yx (solid Vikhlinin ct al. 2008b). All four 
studies are in excellent agreement with each other despite the tight 
cosmological constraints and the different sources of systematic un- 
certainty among the various analyses. 

complement CMB, supernova, and BAO studies. 

We begin by focusing on the somewhat surprising re- 
sult by Vikhlinin ct al. (2008b) that a joint WMAP5 and 
low redshift cluster abundance experiment does not pro- 
duce an interesting constraint on the equation of state 
of dark energy w. The reason this is surprising is that 
WMAP5 has measured the amplitude of the power spec- 
trum at recombination to high accuracy. Given this value 
and a cosmological model, one can predict the value of 
(is today. By demanding that this prediction agrees with 
the cluster normalization condition, one ought to obtain 
a tight constraint on the dark energy equation of state. 

To understand why this is not the case, consider first 
the WMAP5 results. The parameters w and fi TO are 
strongly degenerate given the WMAP5 data alone, as 
shown in Figure 14. The value of erg implied by the 
WMAP5 data depends sensitively on these two param- 
eters, so a large uncertainty in w and Q m dramati- 
cally increases the area of the as — Q m plane allowed 
by the WMAP5 data. Moreover, we can see from Fig- 
ure 14 that the WMAP constraint goes from being or- 
thogonal to the cluster normalization condition to being 
parallel to it, implying that the cluster normalization 
condition cannot improve upon the dark energy con- 
straints of WMAP alone. Indeed, a prior of the form 
cj 8 (fi m /0.25) - 41 = 0.832 ± 0.033 has a minimal impact 
on the error bar in w. Fortunately, given this understand- 
ing, it is easy to see how to improve this situation: we 
need to introduce an additional observable which breaks 
the w — D, m degeneracy. As an example, in the above Fig- 
ures we also show the 68% confidence intervals obtained 
for three additional analyses: 

1. A joint WMAP5+BAO analysis, which includes 
the dark energy constraints derived by Eiscnstcin 
ct al. (2005) using the Baryon Acoustic Oscillation 
(BAO) measurement from the SDSS LRG galaxy 
sample. 
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TABLE 5 

COSMOLOGICAL CONSTRAINTS FROM MULTIPLE CLUSTER ABUNDANCE EXPERIMENTS 



Source 


Reference 


7 


(T g (fW0.25)T 


C8 




maxBCG Richness Function 
X-ray Luminosity Function 
Temperature Function 
Mass function estimated with Yx 


this work 
Mantz et al. (2008) 
Henry et al. (2008) 
Vikhlinin et al. (2008b) 


0.41 
0.62 
0.30 
0.47 


0.832 ± 0.033 
0.85 ± 0.07 
0.80 ±0.04 

0.808 ± 0.024 


0.807 ±0.020 
0.809 ± 0.028 
0.795 ± 0.023 
0.798 ±0.017 


0.270 ±0.019 
0.272 ± 0.026 
0.258 ± 0.025 
0.260 ±0.014 



Note. — The ag and f2 m constraint from the maxBCG + WMAP5 analysis quoted here differs very slightly from 
that presented in Figure 5 because of the simplified approach we have taken in this section for deriving the constraints 
(see text for details). Note Vikhlinin et al. (2008b) quote their result as o-g(n m /0.25)°- 47 = 0.813 ± 0.013 (stat) ± 
0.02 (sys). For this study, we have simply added these two uncertainties in quadrature. 
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Fig. 14. — Parameter constraints on the w — f2 m plane (left) and ag — Q. m plane (right) in a flat u>CDM cosmology, for various data 
combinations. All contours shown are 68% confidence, and are obtained using the MCMC chain outputs downloaded from the LAMBDA 
website (http:/ /lambda. gsfc.nasa.gov/). Despite the fact that the WMAP5 data constrain the amplitude of the primordial power spectrum 
with comparable accuracy in both a ACDM and «>CDM cosmology, allowing w to vary introduces a large degeneracy between w and Q m - 
This degeneracy severely degrades the WMAP constraints in the ag — f! m plane, as seen in the right panel. Adding new observables that 
break the w — tl m degeneracy restores the complementarity between WMAP5 and clusters in the ag — Q m plane, which helps improve dark 
energy constraints through the growth of structure. 



2. A joint WMAP5+SN analysis, which draws on the 
Union combined dataset (Kowalski et al. 2008), a 
compilation of Super Nova (SN) data composed of 
the Ricss et al. (2004, gold sample only), Asticr 
et al. (2006), and Miknaitis et al. (2007) supernova 
samples. 

3. A joint WMAP5+BAO+SN analysis, which adds 
both BAO and SN measurements as extra observ- 
ables. 

In all cases, the confidence contours are estimated 
based on the MCMC data made publicly available 
by the WMAP team through the LAMBDA website 
(http://lambda.gsfc.nasa.gov/). These data sets break 
the w — Q TO degeneracy from the WMAP5 data alone, 
and they restore the complementarity between WMAP5 
and the cluster normalization condition in the og — Q m 
plane, as illustrated in the lower panel of Figure 14. 

Figure 15 shows the constraints in the w — 
o"8(^m/0.25) 41 plane for the various analyses consid- 
ered in Figure 14. The corresponding cosmological con- 
straints are summarized in Table 6. In order to com- 
pute how cluster abundances improve cosmological con- 
straints, we have simply added the cluster normalization 
condition derived within the standard ACDM cosmologi- 
cal model as a prior. While one might worry that letting 
the dark energy equation of state vary would degrade 
the uncertainty in the cluster normalization condition, 
in practice w has a minimal impact provided one re- 
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Fig. 15. — Confidence contours in the to — ag (n m /0.25) 041 plane 
for analyses using various combinations of cosmological data. Con- 
tours all indicate 68% confidence. Provided the w — il m degener- 
acy from Figure 14 is broken by an additional observable, cluster 
abundances can help constrain dark energy through the growth of 
structure between the time of last scattering and the low-redshift 
universe. 

stricts oneself to low rcdshift cluster samples. For in- 
stance, varying w by Aw = 0.1 changes the comoving 
distance to the median redshift of maxBCG clusters by 
w 1%. The growth function is even less sensitive, vary- 
ing by a mere ~ 0.3%. Thus, the cluster normalization 
condition from low-redshift cluster samples is essentially 
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TABLE 6 

Cosmological Constraints in a Flat ujCDM Cosmology 



Experiment Q 

WMAP5 0.266 ±0.086 -1.05 ± 0.34 0.811 ±0.121 

WMAP5±BAO 0.251 ±0.027 -1.20 ±0.24 0.885 ± 0.094 

WMAP5±SN 0.274 ±0.023 -0.98 ± 0.07 0.798 ±0.053 

WMAP5±maxBCG 0.265 ± 0.048 -1.07 ±0.34 0.815 ±0.061 

WMAP5±SN+BAO 0.274 ± 0.015 -0.995 ± 0.067 0.808 ± 0.047 

WMAP5±SN+maxBCG 0.274 ± 0.016 -0.978 ± 0.053 0.801 ± 0.026 

WMAP5±BAO±maxBCG 0.258 ± 0.023 -1.097 ±0.160 0.831 ±0.044 

WMAP5±BAO±SN+maxBCG 0.272 ± 0.013 -0.989 ± 0.053 0.805 ± 0.026 



Note. — The constraints quoted here are derived by multiplying the WMAP5 
likelihoods with a Gaussian prior of our cluster normalization condition. We demon- 
strate in the text that results derived in this way are nearly identical to those from 
more detailed treatments. 



independent of w. 

To demonstrate this explicitly, we compare the re- 
sults of Vikhlinin et al. (2008b) to the results from our 
simple analysis in which the cluster normalization con- 
dition o- 8 (r2 m /0.25) - 47 = 0.808 ± 0.024 from Vikhlinin 
et al. (2008b) is added to the data sets mentioned above. 
For a joint WMAP±BAO±clusters analysis, they find 
w = — 0.97±0.12 while our simple analysis results in w = 
— 1.04 ± 0.13. Note that the error bars are nearly identi- 
cal, while the central value for w differ by only half a a. 
This is still true after adding supernovae as an additional 
constraint, in which case they find w = —0.991 ± 0.045 
compared to our w = —0.971 ± 0.048. A similar conclu- 
sion can be reached for the Mantz et al. (2008) anal- 
ysis. For a joint WMAP±supernovae +/ gas +clustcr 
abundance analysis, the find w = —1.02 ± 0.06. With 
our simple analysis, and ignoring f gas , we obtain w = 
—0.98 ± 0.06. This demonstrates that, to high accuracy, 
current cluster catalogs improve cosmological constraints 
on dark energy only through the low-redshift cluster nor- 
malization condition. 

In summary, we have shown that cluster abundances 
help constrain the dark energy equation of state princi- 
pally through the cluster normalization condition at low 
rcdshifts, which constrains the growth of structure be- 
tween the epoch of recombination and today. However, 
the intrinsic degeneracy between w and f2 m given CMB 
data renders this test ineffective unless the degeneracy is 
broken by an additional data set. It is also worth remark- 
ing here that if we compare the results of a WMAP5±SN 
analysis to those obtained after including the maxBCG 
cluster normalization condition, the constraint on the 
dark energy equation of state is only improved at the 25% 
level, going from Aw = 0.07 to Aw = 0.054. This reflects 
the fact that distance-redshift relationships tend to be 
more sensitive to w than the growth of structure. Nev- 
ertheless, the good agreement between the WMAP5+SN 
constraints and the cluster normalization is far from triv- 
ial. Indeed, the WMAP5+SN likelihood contours in the 
as — &m plane assume general relativity, so the good 
agreement with our data indicates that we are not able 
to resolve any departures from Einstein's theory of grav- 
ity (see also Rapetti et al. 2008; Mortonson et al. 2009). 
While quantitative constraints on such deviations are 
model dependent, our final error on as allow us to unam- 
biguously state that models for which the growth factor 
between last scattering and today differ from our best fit 
ACDM model by « 6% can be ruled out at the 2cr level. 



6.3. Prospects For Improvement 

It is worth considering to what extent we can ex- 
pect the cosmological constraints from maxBCG to im- 
prove with further study. Given that the two principal 
sources of systematic uncertainty are the amplitude of 
the weak lensing mass calibration and the prior on the 
scatter of the mass-richness relation, we focus here on 
those two quantities. More specifically, we re-analyze 
our data using artificially tight priors on each of these 
parameters individually, as well as on both parameters 
simultaneously. The tight priors adopted for this ex- 
ercise are (3 = 1.00 ± 0.01 for the mass bias parame- 
ter, and ctj\/|at 200 = 0.45 ± 0.02 for the scatter in mass 
at fixed richness. We find that the uncertainties in 
5*8 = (78(^m/0.25) a41 for each of these analyses are: 

1. fiducial analysis: ASs = 0.033, 

2. tight scatter prior alone: ASs — 0.026, 

3. tight mass bias prior alone: ASs — 0.025, and 

4. tight mass bias and scatter prior: ASs = 0.018. 

Thus, a tight prior on either the mass bias or the scatter 
parameter improves the principal maxBCG cosmologi- 
cal constraint by 25%. If both priors are tightened, the 
improvement is as high as 50%. The corresponding con- 
straint on w for this most optimistic scenario, assuming a 
joint WMAP5+BAO±SN+maxBCG analysis, would be 
Aw = 0.049, which is only a 7% improvement relative to 
the current constraint. 

Are such improvements feasible? In principle, yes. Im- 
provement of the mass bias parameter is possible through 
a follow-up spectroscopic program aimed at calibrating 
the mean lensing critical surface density of the lens- 
source pairs used to estimate the mean cluster masses. 
Likewise, an extensive X-ray follow-up program could in 
principle constrain the scatter in the mass-richness re- 
lation to high accuracy. In practice, realizing such tight 
priors might be difficult. For instance, given the current 
scatter estimate &m\n 20 o = 0-45, we require « 400 X-ray 
follow ups to achieve an uncertainty of AaM\N wa ~ 0.02. 
Such an extensive program seems unlikely to be feasi- 
ble any time in the near future. What is needed, then 
is a way to significantly reduce the number of follow- 
up observations necessary to improve our cosmological 
constraints. Fortunately, the number of follow-up obser- 
vations necessary to achieve such accuracy scales as the 
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square of the scatter <Jm\n 200 , so the best thing to do at 
this point is probably to focus on constructing new rich- 
ness estimators that better correlate with halo mass (see 
e.g. Rozo et al. 2008b; Reyes et al. 2008). 

Finally, we note that there remains additional informa- 
tion about the maxBCG clusters that has not yet been 
incorporated into our analysis. This includes galaxy ve- 
locity dispersion data (Becker et al. 2007), X-ray data 
(Rykoff et al. 2008, while X-ray data was used to place 
a constraint on the scatter in mass at fixed richness, we 
did not otherwise use the X-ray data in this analysis), 
and clustering information (Estrada et al. 2008). Includ- 
ing these additional probes of cluster mass should help 
further improve our cosmological constraints. 

7. SUMMARY 

We have performed a joint analysis of the abundance 
and weak lensing mass estimates of the maxBCG clusters 
detected using SDSS imaging data. In addition to this 
data, a prior on the scatter in the mass-richness relation 
derived from demanding consistency between the weak 
lensing and X-ray mass estimates of the clusters. Our 
cosmological constraints can be summarized as 

cr 8 (fi m /0.25) - 41 = 0.832 ± 0.033, (18) 

which is consistent with and complementary to the latest 
WMAP results. With a joint maxBCG and WMAP5 
analysis we find 

cr 8 = 0.807 ±0.020 (19) 
n m = 0.265 ±0.016. 

These results firmly establish optical cluster studies as a 
method for deriving precise cosmological constraints. Im- 
portantly, our results are in excellent agreement with and 
of comparable precision to X-ray derived cluster abun- 
dance constraints, clearly demonstrating the robustness 
of galaxy clusters as a tool of precision cosmology. 

We have discussed how and why galaxy clusters can 
help constrain dark energy evolution, demonstrating that 
even in those data sets where the evolution of cluster 
abundance with redshift is clearly detected, constraints 
on the dark energy density and equation of state of a 
joint WMAP and cluster abundance analysis are domi- 
nated by the low redshift cluster normalization condition 
crgfi^ = constant. These joint constraints are driven by 
the growth of the matter fluctuations between the time of 
last scattering and the low redshift universe. Thus, while 
cluster abundances provide only moderate improvements 
to dark energy constraints derived from joint WMAP and 
supcrnovae analysis, we have argued that they provide 
an important consistency test of general relativity. More 
specifically, our constraint on eg allows us to rule out at 
the 2cr lavel any models for which the growth of struc- 
ture between last scattering and today differs from that 
of our best fit ACDM model by more than ~ 6%. 



At this time, the dominant systematic uncertainty in 
our analysis is the uncertainty in the weak lensing mass 
scale due to scatter in photometric redshift estimates of 
the source galaxies. In addition, improvements to our 
understanding of the scatter of the mass-richness relation 
could help tighten our cosmological constraints. Follow- 
up observations can help in this regard, but the number 
of follow ups necessary to have a significant impact on our 
results is currently very large. Fortunately, reducing the 
scatter of the mass-richness relation (see e.g. Reyes et al. 
2008; Rozo et al. 2008b) may help reduce the number 
of follow up observations necessary to achieve improved 
constraints. 
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